IQSL-Grating Cavity 



Strong Optical Confinement between Non-periodic Flat Dielectric Gratings 
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We present a novel design of optical micro-cavity where the optical energy resides primarily in 
free space, therefore is readily accessible to foreign objects such as atoms, molecules, mechanical 
resonators, etc. We describe the physics of these resonators, and propose a design method based on 
stochastic optimization. Cavity designs with diffraction- limited mode volumes and quality factors 
in the range of 10^-10^ are presented. With a purely planar geometry, the cavity can be easily 
integrated on-chip using conventional micro- and nano- fabrication processes. 

PACS numbers: 42.60.Da, 42.15.Dp, 42.25.Fx, 87.80.Cc 



Optical cavity has become an indispensable element 
in a wide range of appHcations such as lasers, sen- 
sors, cavity quantum electrodynamics (CQED), op- 
tomechanics, optical traping and manipulation, etc.[l|- 
0]. For these applications, a high quality factor (Q) 
and/or a small effective mode volume (defined as F = 
J e\E\^dY/ {e\E\^/2)^^J[8] is usually important, with 
the first a measure of the photon temporal confine- 
ment, and the second inversely proportional to the elec- 
tric field per confined photon. The ratio of the two 
(Q/V), which is proportional to the Purcell factor [oj, 
represents the increase of the spontaneous emission rate 
for a light emitter and is another figure of merit in cav- 
ity design. It plays an important role in determining 
the lifetime of an excited atom or a quantum dot con- 
tained within the cavity. Various on-chip optical cavi- 
ties of micrometer size have been designed and fabricated 
that can be categorized into three types, depending on 
their light-confinement mechanisms: Fabry-Perot (FP) 
cavities, cavities with whispering-gallery (WG) modes 
(micro-rings, disks, spheres, toroids, etc.) and two di- 
mensional (2D) photonic crystal (PhC) defect cavities [l|. 
Integrated FP cavities — such as those used in Vertical 
Cavity Surface Emission Lasers (VCSELs) — go to reso- 
nance when light interferes constructively after a round 
trip between two mirrors that are usually distributed 
Bragg reflectors (DBRs). High Q-s close to 10^ are 
reported, but a current-confining oxide aperture or a 
waveguide structure such as a micro pilla r is usually re- 
quired for lateral mode confinement |1Q[ [ii| . while the 
mode often penetrates into the DBRs for a depth up to 
several wavelengths. WG cavities confine the light along 
a circular or spherical wall through total internal reflec- 
tion (TIR), and build up a resonant mode when the phase 
delay in a full circular trip is the multiple of 27r. Very 
high Q values (up to 10^) were achieved for devices of 
large size 100/im) with less ideal V^[l]. Micro-rings 
of small radius were also reported with the trade off of 
lower Q-s due to a larger bending loss [12]. 2D PhC de- 
fect cavities confine the light through a photonic bandgap 
in the lateral directions and TIR in the vertical. They 
are capable of achieving modes with small mode volumes 



1 1 1 inimF' ■ ■ ■ ■ 1 1 




FIG. 1. (a) Plot of the instantaneous distribution of Ex (hor- 
izontal in the plane of the page) reflected by a parabolic 
grating mirror, normalized to its maximum value. The in- 
cident field has a flat curvature and an operating wavelength 
Ao = 835 nm. (b) Plot of Ex at 835 nm for a cavity composed 
of two such mirrors, normalized to its maximum value. The 
white scale bars have a length of 1 /xm. 



(~ (Ao/n) ) but usually lower Qs compared to the best 
of WG cavities, mainly due to the difficulty in confining 
the mode in the vertical direction when the size is small 



One essential issue in many applications using opti- 
cal cavities is the accessibility for particles to regions of 
high electric field. This determines the coupling strength 
between the light emitter and the cavity in CQED, the 
sensitivity of an optical-cavity based sensor, and the 
strength of the force when the cavity mode is used to 
influence nano-mechanical structures or to trap and ma- 
nipulate nano particles. While each of the cavity designs 
discussed above manage to achieve high Q and/or small 
in reality they all show disadvantage in this crite- 
rion. For WG and PhC cavities that utilize TIR, the 
high-field regions are either contained inside the dielec- 
tric or in a small region close to the dielectric sidewalls, 
making the access to this region very difficult. In reality, 
the atoms, molecules or nanomechanical resonators are 
either brought into the evanescent tail of the mode out- 
side of the cavity [it 0, [gI, or planted inside the dielectric 
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body during the fabrication [H^. Although an FP cavity 
without any medium in between the two mirrors could 
in principle provide a mode accessible to foreign parti- 
cles, the lack of lateral confinement in this case would 
inevitably result in a largely extended, or even instable, 
resonant mode[l^. 

In this paper we present a new type of optical cav- 
ity with a reasonably high Q, a diffraction-limited mode 
volume, and much-improved field accessibility. To intro- 
duce our cavity design, let us first examine a conventional 
FP cavity, but consists of two identical dielectric mirrors 
working as parabolic reflectors. Following an approach 
we recently developed in Ref. Il5|, these dielectric mir- 
rors were designed by introducing a smooth variation of 
the duty cycle (DC, the ratio of the width of the dielec- 
tric groove to the period of the grating) of a dielectric 
grating. By locally tuning the DC, the reflection phase 
of the mirror is tailored to be a parabolic distribution 
to the position while the reflection magnitude maintains 
high everywhere, so that the mirror works as a parabolic 
reflector while remaining a planar geometry [iBj. Fig. [Tal 
shows a 2D finite difference time domain (FDTD) simu- 
lation of the behavior of such a mirror under a Gaussian 
beam incidence. This simulation (and others below) was 
done using an open source package MEEP [l6| . The mir- 
ror has a focus length of 13/im and a reflectance of 98%. 
The resonant mode of a cavity composed of two of such 
mirrors is shown in Fig. [lb] together with the outline of 
the geometry. Notice that, the mode has a Q = 2,000 
with a full width at half maximum (FWHM) of I.SAq. 
These can not be achieved for VCSEL cavities of planar 
reflectors without a mode- confining mechanism. 

Compared to conventional VCSEL cavities with DBRs, 
the one shown in Fig. [lb] supports a strictly stable mode 
of standing Gaussian beam with a finite cross section, 
with a ~ l/4Ao field penetration into the dielectric mir- 
rors. It is in principle similar to a conventional FP res- 
onator with curved mirrors but is compatible to a pla- 
nar fabrication process. The relevant parameters such as 
the resonance frequency, the waist, and the quality fac- 
tor can be predicted from the effective curvature of the 
mirrors and their separation using formulas for FP with 
curved mirrors 0. However, this approach would not 
be successful when we try to further reduce the mode 
cross-section to sub-wavelength size. A Gaussian beam 
with a cross section even smaller would be too diverging 
to satisfy the longitudinally slow- varying condition, thus 
the underlying physics of the cavity shown in Fig. [lb] no 
longer holds. Also, the optical mode may only interact 
with a few individual dielectric "grooves" (or "posts"), 
whose collective behavior as a grating is no longer rel- 
evant. Rather, the dielectric grooves would behave like 
independent, near-resonant scatterers, and the global be- 
havior of the assembly as a cavity is the result of the com- 
bined effect of these coupled scatterers. If the structure 
could be designed so that the radiation from the scat- 
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FIG. 2. (a) The schematic of the cavity design, (b)-(d) The 
instantaneous field distributions of one cavity design, normal- 
ized to their respective maxima: (b) 1^1^; (c) Hz; (d) Ex. The 
white scale bars in these plots have a length of 1 /im. (e) \E\'^ 
distribution along the horizontal and vertical symmetry axis. 



terers interferes destructively in the far-field, a high-Q 
optical mode can be obtained. 

A deterministic design of the resonance features of the 
grooves and their mutual coupling would be extremely 
complicated. As a first attempt to test the feasibility 
of the idea we use a stochastic optimization method to 
search for a cavity design with high Q and/or small V. 
We first try to design an optical cavity with a structure 
symmetric with respect to the x and y axes (Fig.[2a|). and 
we look for a fundamental TM mode near Aq = 1.55/im 
with an asymmetric transverse magnetic field distribu- 
tion along the x axis. In each quadrant, we use N dielec- 
tric grooves of refractive index n = 3.48 (silicon, around 
Ao = 1.55/im), where N is determined before hand. The 
width and position of each groove is determined by pa- 
rameters ai and bi (as shown in Fig l2aj) . Some or all of the 
ai and 6^, together with the groove thickness t and cavity 
length L are used as optimization parameters. This cor- 
responds to 14 degrees of freedom for A/" = 6. An exhaus- 
tive search is not practical, and a stochastic optimization 
has proved to be essential here. A parallel genetic algo- 
rithm (GA) is used in our study, and we refer the readers 



to Ref. [ITI for a brief, yet informative introduction to 
GA and its applications to EM problems. In a GA pro- 
cess, each optimization parameter is binary-coded and 
the codes are then concatenated to form a binary string 
(a "chromosome"). Each chromosome represents a possi- 
ble cavity design whose behavior is simulated by an EM 
solver and its "fitness" , a quantitative description of how 
well the design approaches the desired performance, is 
calculated. The GA evolves towards a design with max- 
imum fitness. Different choices of the fitness function 
give preference to different aspects of performance, thus 
resulting in different cavity designs. In our study we use 
either Q or Q/y as the fitness function. The EM simula- 
tion is done on MEEP employing its parallel algorithm. 
The optimization process is carried out on a computer 
cluster using up to 160 cores, with multiple cavity simu- 
lations, each of which uses 4 cores, running at the same 
time. Because of this massively parallel strategy, the op- 
timization process converges quickly ('-«^20-30 hours) in 
spite of the large search space. 

The optimization process indeed yields cavities with 
properties as expected. The normalized instantaneous 
distributions of i^^, and Ex for one of the results 

are shown in Fig. [2b| [2c| and [2dl respectively. In this 
case, the thickness of the grooves is t = 400 nm and 
the cavity length is L = 759nm. The structure has a 
total width of about 8/im and supports a cavity mode 
with a free space wavelength Aq = 1542.6nm with a qual- 
ity factor Q = 14, 697. Notice that the profile of the 
grating-like structure is fairly "irregular" compared to 
that of the parabolic mirror shown in Fig. [TJ The cav- 
ity does show reasonable tolerance of variations in the 
structural parameters around the optimized design. For 
example, Q decreases as L deviates from the optimal 
value, but stays above 10, 000 for 750nm < L < 770nm, 
when the resonant wavelength varies between 1539nm 
and 1548nm. The cavity mode, which otherwise would 
not exist, demonstrates the radiation cancellation pro- 
vided by each dielectric groove, and the quality factor is 
determined by the residue radiation into the far field that 
is not ideally canceled by the grooves. 

In addition to the fact that it can be built using planar 
technology, two other features of the cavity deserve spe- 
cific mentioning. First, the mode is concentrated around 
the center of the structure confined in free space between 
the dielectric layers. As shown in Figl2el along the y axis 
the field is confined primarily between the two dielectric 
gratings with a full width at half maximum (FWHM) 
of only ~ 400 nm. Along the x axis, where no mode- 
confining structure exists, the FWHM is ~900 nm, well 
below the vacuum wavelength. The equivalent volume 
for the 2D mode is V = 0.283 /im^ = 0.119 Ag. Sec- 
ond, the radiation leaked into the far-zone forms multi- 
ple sharp beams, with the dominant beam pointing to 
±7/ directions for many cases. For the radiation, the di- 
electric grooves work as an antenna array of large aper- 
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FIG. 3. (a) A schematic of the 3D cylindrical cavity, (b) The 
instantaneous distribution of Er (top), Ez (center) and \E\'^ 
(bottom) of the cavity mode, all normalized to their respective 
maxima. The white scale bars have a length of 1 /xm. 



ture fed by the cavity mode of small cross section. By 
reciprocity, we can also use this "antenna array" at the 
receiving mode, and to excite the cavity mode using a 
plane incident wave. This is indeed proved by numerical 
simulation. This type of excitation is much more conve- 
nient compared to the excitation of WG cavities where 
waveguides and phase-matching are required, or that of 
PhC cavities where tapered fibers are used. 

The simulation cases shown so far are for planar, 2D 
structures invariant along one of the Cartesian axes. We 
can also design 3D cavity that has a cylindrical sym- 
metry. In this case the cavity consists of two identical 
layers comprising multiple concentric rings of different 
widths, as schematically shown in Fig. [3al and we look 
for solutions in which the angular dependence of the fields 
is proportional to e*^"^, where (j) is the azimuthal angle 
and m is an integer. Like the linear grooves, cylindrical 
grooves exhibit local resonances that can be harnessed 
to support localized 3D cavity modes. Thanks to the 
rotational symmetry of the structure the FDTD simu- 
lation reduces to an effective 2D problem that can be 
quickly solved. We apply the same GA optimization to 
the cylindrical problem using the thickness, widths and 
radius of the dielectric rings as optimization parameters. 
In what follows we limit our analysis to TM modes with 
m = (i.e. the only non-zero component of magnetic 
field is H(f)). Cavities supporting m ^ modes can also 
be designed using a similar approach. Fig. [3b] shows in- 
stantaneous distribution of E^^ E^^ and on the cross 



4 



□ 


□ □ □•□a □ □ 


□ 














□ □ □ □ 


□ 






(a) 




(b) 1 


D 




D 


□□ m □□ 1 


D 


I □□man d: 


D 




(c) 


(d) 



FIG. 4. \E\'^ distributions (normalized to their respective 
maxima) for various cavity designs, (a) 2D cavity, limited 
size, regular groove shape, Q = 78, 260. (b) 3D cylindrical 
cavity, higher order mode, Q = 92, 384. (c) 2D cavity, thick 
dielectric layer, Q = 601, 320. (d) 3D cylindrical cavity, open 
central part, Q = 14, 337. The white scale bars have a length 
of 1 /xm. 



section through the symmetry axis for a typical design of 
Ao = 1552.6 nm and Q = 60,435. Similar to the 2D cav- 
ity case, the mode is largely confined in free space, around 
the symmetry axis in between the dielectric layers, with 
an effective mode volume of V = 0.637 /im = 0.264 Xl 
The cavity can also be excited through an azimuthally 
polarized beam at normal incidence. 

The examples shown above are only a subset of the de- 
signs we have achieved. With different choices of fitness 
function (Q, l/V, or Q/V), types and ranges of the op- 
timization parameters (groove thickness, cavity length, 
some or all the sizes and positions of the grooves, etc.), 
cavities with widely varying characteristics can be de- 
signed. A few other typical examples are shown in Fig. 31 
where Fig. 2^ and c show 2D designs while those in 
Fig.|4)3 and d are cylindrical 3D cavities. In the design of 
Fig. 2^ we limit the over all lateral size of the cavity to be 
smaller than 6.5 /im and use regularly shaped dielectric 
grooves with aspect ratio < 1, to reduce the fabrication 
difficulty. The mode is wider compared to that in Figl2] 
but the effective 2D mode volume is still smaller than Aq 
with Q > 78, 000. For the design in Fig. St, thicker dielec- 
tric grooves are used and more EM energy is stored inside 
the dielectric, resulting in a quality factor Q > 600, 000. 
In Fig. HJd a cavity of a higher order longitudinal mode 
with a quality factor Q > 90, 000 is shown. In the cavity 
in Fig. |4]i a large portion in the center of the structure 
is left open for convenient access of the mode field. The 
optical field is still concentrated around the symmetry 
axis with Q > 14, 000. 

We emphasize here that the function of the two layers 
of dielectric grooves in this cavity design is fundamen- 
tally different to that of the mirrors in a FP cavity. A 



single layer of this structure does not show any appre- 
ciable reflectivity, and the field in between the two lay- 
ers can not be described as a back-and-forth traveling 
wave. As we see in the plots, a large portion of the res- 
onant mode, including the maximum of the field, stands 
in free space, freely accessible by molecules, light emit- 
ters, nanomechanical resonators, etc. This would provide 
a much stronger effect compared to the evanescent field 
coupling used in existing references. The cavity holds a 
resonant mode with an open structure, essentially works 
as a ca ge" for photons. Because of these features, we 
believe it can be a unique candidate for a wide range of 
applications, many of which are currently under study. 
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